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TECHMCAL NOTE 5O77 


THE EFFECT OF DYNAMIC LOADING ON THE STRENGTH 
OF AN INELASTIC COLUMN 

By William A. Brooks, Jr., and Thomas W. Wilder, III 


SUMMARY 


The maximum loads of idealized inelastic H-section columns whose 
pinned ends approach eaoh other at a constant rate are presented. The 
solutions indicate that as the rate of end displacement heccmes smaller 
the dynamic buckling solutions approach the static solution as a lower 
limit. The effects of inertia forces are apparently negligible at rates 
of end displacement comparable to those normally used in static column 
tests. For all rates of end displacement investigated the static maxi- 
mum load may he employed as a conservative estimate of the maxi mi mi 
colimm load. 


INTRODUCTION 


Only in recent years has the solution to the problem of the stat- 
ically loaded Inelastic column been siifficlently developed to provide 
answers as satisfactory as those predicted from elastic-column theory. 
Shanley (ref. l) gave new life to inelastic -colimm theory by Indicating 
that a column may start to bend at the tangent -modulus load and that 
the maximum load of a straight inelastic column is always less than the 
reduced-modulus load but greater than the tangent -modulus load. By 
extending Shanley 's theory, Duberg and Wilder (ref. 2 ) found that the 
maximum load could be determined more definitely and that it depended, 
for a large part, on the shape of the stress -strain curve. It was later 
shown (ref. 5) that the maximum load of an Initially curved inelastic 
columin may be less than the tangent-modulus load, and that initial out- 
of -straightness has a greater effect on the maximum load of an inelastic 
column than on that of an elastic column. 

V/hereas the behavior of columns subjected to static loading condi- 
tions is comparatively well known, investigations involving the behavior 
of dynamically loaded col\mms have been limited in scope and mmiber in 
spite of the fact that most structinral loads, in particular the most 
significant aircraft structirral loads, are dynamic. Hoff and his asso- 
ciates (refs. ^ and 5) treated the case of an elastic column whose ends 
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approach each other at a constant rate and showed that for rates of end 
displacement ccm 5 )arahle to those obtained in static-testing machines 
there is good agreement between the maximum load and the Euler load. It 
was also shown that at higher rates the maximum load becomes greater than 
the Euler load. The Inelastic colimm was analyzed in much the same man- 
ner by Chawla (ref. 6), who foimd that for relatively slow rates of end 
displacement the Tna-x-iTmTTn column load was much smaller than the tangent- 
modulus load. Chawla 's results appeared low; therefore, a static solu- 
tion was made by the authors which revealed that the maximum load obtained 
in reference 6 was not only less than the tangent-modulus load but also 
less than the static maximum load.. If this siorprislng result is correct, 
the conseqiiences are indeed serious because its indication is that a 
static analysis could not be used to obtain a conservative estimate of 
the maximum load. 

The present paper gives the results of an analysis which was made 
to ascertain whether some particular combination of the factors which 
influence the strength of inelastic columns would cause the strength of 
the column to be less when rapidly loaded than when slowly loaded. The 
only manner of loading which is considered is that which causes the 
pinned ends of the column to be axially displaced with a constant veloc- 
ity. The factors considered most imiportant for the study are the ini- 
tial out-of-strai^tness, the shape of the stress -strain curve, the 
slenderness ratio of the column, and the rate at which the ends of the 
column are displaced. The parameters which describe these factors are 
varied to yield numerous solutions which show the inflixence of each 
factor on the strength of the column. 

The column used in the analysis is the idealized H-sectlon column. 

The analytical expression used for the stress -strain curves is that 
proposed by Ramberg and Osgood (ref. j) and is assumed to be unaffected 
by the rate of loading. 


SYMBOLS 


A cross-sectional area of column 

b colxmm thlcJmess (see fig. l) 

d total lateral deflection of column at midheight after 

loading 

initial lateral deflection of colimm at midhel^t 
before loading 



MCA TN 5077 ■ 


3 


e 


E 

Eg 

f 

L 

m 

n 


q. 

t 

V 

X 

y 

yo 

€ 

^1 


^E 


I. 

p 

a 

^1 

0 "^ 


dimensionless initial lateral deflection of column at 
midlifiight before loading, 2dQ/b 

Young's modulus 

secant moduliis 

dimensionless total lateral deflection of coltrain at 
mi dhe l gbt after loading, 2d/b 

column length 

mass per unit volume of column material 

Ramberg-Osgood stress-strain-curve shape parameter 
(see ref. 7) 

lateral inertia force per unit length 
time 

velocity at which the ends of the column are displaced 
longitudinal distance measured from end of column 
total lateral deflection of column 
initial lateral deflection of coltmm 

strain 

elastic strain corresponding to stress 0^ 

elastic strain corresponding to Euler stress, 

dummy variable of integration 
radius of gyration, b/2 
stress 

O.7E secant yield stress 

cross-sectional stress corresponding to tangent- 
modiilus load 
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a average cross-sectional stress of column 

T dimensionless time parameter^ vt/Le^^ 

CD first natural circular frequency in transverse bending 


n 


dynamic buckling index. 



Subscripts : 

L left (concave) flange 

R right (convex) flange 


max Tnn.ylTTTiTm 

rev corresponding to the strain at which reversal occurs 

Dots indicate differentiation with respect to time t; primes 
indicate differentiation with respect to dimensionless time t. 


THEORY 


The idealized pin-ended colimm (fig. 1) treated in the present 
paper consists of two thin flanges of equal area separated by a web of 
Infinite shear stiffness and of negligible area. Inasmuch as the veloc- 
ity of compression-wave propagation is large compared with bcth the rate 
of end displacement and the rate of the resulting lateral deflection, a 
state of uniform average compression is assumed to exist along the length 
of the column. 

The Initial conditions at time zero are 




y = Yo 

y = 0 
V = 0 


j 


( 1 ) 


If D'Alembert’s principle is applied to the portion of the colimon 
shown in figinre 1, the following equilibrium equation can be obtained 
by summing the moments about the pin: 
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This equation relates the flange stresses and to the lateral 

deflection y. Another relation between the stresses and latereil deflec- 
tion can be obtained through the use of the stress -strain relation of 
the material In conjunction with the satisfaction of ccm^atlblllty between 
the flange strains and the lateral deflection. 

In the present paper the stress-strain curve Is assimied to be rep- 
resented by the Ramberg-Osgood equation: 


_L = _2_ + 5 , 
^1 ^1 7 (^ 1 ] ^ 5 ) 

This equation Is used when the material Is loading. Unloading (strain 
reversal) which may occur In the convex flange Is assumed to take place 
elastlcallyj thus 


€ _ cr 5 /,^;e^ 

€i 01 7 \ tJi / 


(if) 


where 0 rev stress corresponding to the strain at which reversal 

occurs . 

Satisfaction of ccmipatlblllty requires that 




b - yp) 

2 dx^ 




6 


. t d^(y - yp) 

2 




(5) 


where, by virtue of the assumption of uniform average compression, the 
average strain "g Is related to the end displacement as follows; 
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Le = vt + 





dx 


( 6 ) 


By assuming the initial out -of -straightness to he 


yo = sin ^ (7) 

and the total deflection to he 

y = d sin 3 ^ ( 8 ) 

1j 


and hy satisfying equilibrium only at midhei^t (x = L/ 2 ), equation ( 2 ) 
becomes 




(I - - ^(1 + a) + ^ ' 0 


( 9 ) 


IThen equations ( 6 ), (7)^ and ( 8 ) are \ised^ the strains (5) become 


'-B. 




( 10 ) 
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Employment of the relations 


- 


it- 2 
^ L 


f = 




e = 


h 


2 


reduces equations (9) and (lO) to 




(u) 


and 


= L 


^ - i Cj.(f2 - e2) H. cj.(f - e) 


6 - 
- IT 


^ - t "Elf" " ^ 


(f 2 - e^) - ej.(f - e) 


( 12 ) 


The differential equation (U) can he made (dimensionless by using 
the relation 


(T^ = E€3_ 


(13) 


and the parameter 


T = 


vt 

Le. 





8 


KACA TN 5077 


A value of t equal to xmlty corresponds to an. end displacement suf- 
ficient to produce a strain of € 3 _ in a straight coltmin. Equation (U) 
thus hecomes 


f - 2t!_ ^ 
”2v2 


^ (f - 1) + ^ (f + 1) 




^1 


( 15 ) 


Equations (12) axe made dimensionless by dividing both sides of 
the equations by and \ising equation (l4) : 




(16a) 


(16b) 


The term before the brackets of equation ( 15 ) uiay be written as 


2 v^ 



( 17 ) 


where, by definition. 


a = 



(18) 


In the expression for ft, the quantity o) is the first natural circtilar 
frequency in transverse bending. The differential equation ( 15 ) may now 
be rewritten: 


f" = 


ft 


-(^e/^i)' 


^(f - 1) + ^(f + 1) 


( 19 ) 
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Simultaneous solution of equations (5) (or (3) and (^ 1 -)), (16), 
and (19) yields the desired lateral deflection. 

An additional equation which proves useful in a numerical solution 
is obtained by using the relation 


_ 2 _ = ^ 

01 E 61 


(20) 


and combining equations (16) and (19) to produce 


f" = 2 _ 


(^e/^i) 


Ec 


E 


T +S(f -a) 

IT 61^ 61' ^ 


(f - 1) + 


E 




T - 1 !S(f 2 . e 2 ) - S(f - e~) 

L ‘*'1 '1 J 

(f + 1 ) ’• 


( 21 ) 


A numerical solution of the problem is described in detail in the 
appendix. 


RESULTS AND DISCUSSION 


The cases for which solutions were obtained are given in table 1 . 
The parameters which are varied are; 

(1) The shape of the stress-strain curve^ which is determined by 
the value of n assigned to the Ramberg -Osgood equation 

(2) The slenderness ratio of the col-umn^ which is determined by 

the value assigned to (which is, in turn, defined by the valvie 

assigned to 0 t/ 0 i) 

(5) The initial out -of -straightness, which is determined from the 
vailue of the dimensionless parameter e 

(il-) The dynamic -buckling index D, which is inversely proportional 
to the square of the velocity at which the ends of the column approach 
each other. 
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Aa explanation of the meaning of the djniamic -buckling index £l can 
he had by realizing that when a paxtictilar column is being loaded, the 
end displacement will reach a value equal to e-gL in a time equal to 

that required for an identical unloaded column to undergo a natural 
vibration of \[fl/2rt cycles. Therefore, 0. is the parameter which can 
be varied to correspond to vario\is rates of end displacement. Table 2 
indicates the approximate range of the velocities employed. 

The dimensionless parameter t (eq. (l4)) may be considered to be 
a shortening parameter, as the product vt is the amount of shortening 
that occvirs in time t. With the exception of figiire 2, the results eire 
presented as plots of load against lateral deflection or shortening. 

Plots of dimensionless lateral deflection against dimensionless 
shortening are given in figure 2 for n = 10 and e = 0.01. The cirrves 
result from the dynamic solutions, and the symbols in this figure rep- 
resent points in the neighborhood of the maxi mi mi load obtained from the 
static solution. At the higher velocities (small values of ft), for a 
given amount of shortening, the lateral deflection is considerably 
smaller than that obtained from the static solution. Th\is, in the 
dynamic solutions, more sho 2 rbening initially resiilts from axial com- 
pression than from bending, and the lateral deflections initially lag 
behind those obtained from the static solution. However, as time pro- 
gresses it is possible for the lateral deflection obtained for all rates 
of end displacement to overshoot that obtained from the static solution. 

For ft = 10^, the lateral deflections from the dynamic solution oscillate 
about those obtained from the static solution. The oscillations are so 
small that it is difficult and impractical to show them. For a.i 1 prac- 
tical purposes, the curve for ft = 10 ^ coincides with that obtained 
from the static solution, and the effects of inertia forces are negligible 

Figure 3 presents the variation of load with lateral deflection for 
the cases given in figure 2. The symbols in this figure and subsequent 
ones represent points in the neighborhood of the maximu m load obtained 
from the static solutions given in reference 5- It is to be noted that 

maximum loads are not indicated for ft = 1 , 10, and 10^. For these 
cases the solutions were terminated when strains became approximately 
equal to 2 percent. Excessive deformation rather than maximum load is 
believed to be the governing criterion for cases in this range of 
velocities . 

In figure h, load is plotted against shortening. The material 
stress-strain curve is also shown and, as would be expected, serves 
as the upper limit. However, it is important to realize that at the 
hlg^r velocities the material stress -strain curve may be altered. In 
the present paper, the rate of end displacement is assumed to have no 
influence on the shape of the stress -strain curve. 
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In order to investigate the effect of initial out-of -straightness, 
a groijp of solutions was made for a given rate of end displacement and 
various values of Initial cirrvature. The solid cirrves of figures 5 
and 6 are these solutions and the dashed curves are the corresponding 
static solutions. As would be expected, the largest values of load 
overshoot are associated with the smallest initial out-of -straightness - 
a situation similar to that which exists for static loading; namely, the 
smallest eccentricities permit the column to bear the largest loads. 
However, it is apparent that initial curvatures, particularly the small 
out-of -straightness, have more effect on the maxiTnum load for dynamic 
loading than for static loading. (The static maximimi load for e = 0 
is practically coincident with that for e = 0 . CXDOOl . ) 

Figure 7 shows the results of solutions for colimins of different 
lengths and of a given material, initial out-of -straightness, and rate 
of end displacement. For coniparison purposes the tangent -miodulus load 
and static maximum load are shown. The static curves were obtained 
from reference 3* A static curve is not given for (or 

equal to 1.1 because, as a result of computational difficulties, this 
case was not considered in reference 3* From these plots it is evident 
that the tangent-mod'ulua load may in some Instances, depending on the 
column proportions, material, and initial out-of -straightness, be an 
unconservative approximation to the maximum dynamic load. For the case 
shown, the tangent-modulus load is a better approximation for the dynamic 
maximum load than for the static maximum load. 

Figures 8 and 9 show solutions for a different material (n = 2). 

A carcparison of these figures with figures 3 ^ reveals that, for a 

given rate of end displacement, the percentage increase in maximum load 
above that predicted by the static analysis is greater for the material 
having a more rapidly cirrving stress -strain curve (n = 10). The load 

overshoot for = 10^ is, in the given case, approximately 32 percent. 
Considering the fact that the initial out-of -strai^tness is rather 
large, it is obvioiis that for n = 2 the dynamic maximum loads for 
the smaller values of initial out-of -straightness may be considerably 
larger than the static maximum load. However, the static maxiirrum load 
is still on the conservative side of the dynamic maximum load. 


CONCLUDING REMARKS 


The solutions presented indicate that as the rate of end displace- 
ment becomes smaller the dynamic buclsling solutions approach the static 
solution as a lower limit. The upper limit for dynamic buckling is, as 
would be expected, the material stress -strain curve, which may be altered 
by the rate of loading. The effects of inertia forces are apparently 
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negligible at rates of end displacement comparable to those no rmall y 
used in static column tests. From the results obtained it can be con- 
cluded that the static maximum load may be employed as a conservative 
estimate of the mp.-xiTTTinn load of a col umn j regardless of the rate of 
end displacement. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., December 29, 1955- 
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APPENDIX 

NUMERICAL SOLUTION 


The nonlineax natiare of the differential equation of lateral column 
motion makes it necessary to resort to a numerical solution. The solu- 
tion is facilitated by the fact that the equation of motion does not 
contai n the velocity of lateral deflection; only the deflection pmd 
acceleration are involved. By approximating the acceleration with an 
interpolation formula and integrating twice, a recursion formula for 
forward integration can he obtained. Once the type of integrating for- 
mula has been decided upon, it is necessary only to obtain the required 
initial values and to choose the interval size. In the sections that 
follow, each of the aforementioned items will be discussed in detail. 


Reciirsion Fomulas for Forward Integration 

When f is known, the strains and ep can be found by sub- 

stituting f and the corresponding value of r into equations (16). 

The stresses Cp and cp can then be determined by unlrig the stress - 
strain relation (3) or relations (3) and ( 4 ). If the stresses that exist 
at f are known, f" can be computed from equation (19). Equation (19) 
may therefore be written in functional notation as 


f" = g(f,T) (Al) 

and can be integrated by \ise of the following recursion formula (ref. 8 ): 



However, as the following diagonal -difference table indicates, 
A^f"p_2_ is a fionction of f and is therefore indeterminate until 
f"p+iy which is, in turn, dependent on fp+p (eq. (Al)), has been 
found. 
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DIAGOML-DUFEERENCE TABLE 


T 

f 

f" 

Af" 

A^f" 

(r - 2 )At 

^r-2 

f " 

^ r-2 

Af" _ 





r-2 


(r - 1 )At 

^r-1 

f " 

^ r-1 

^"r-1 


(r)AT 

fr 

f" 

r 

Af"^ 


(r + 1 )At 

^r+1 

-pii 

^ r+1 




An approximate solution (ref. 8) which does not require 
is the following; 

fr+i = 2fj. - f^.i + (At)2 ^ ^^"r- 2 ) (A5) 

With the approximate value of f^+i given hy equation (A5), an approxi- 
mate value of A2f"^_^ may he found. Equation (A2) is then used to 
correct and check the approximate valvie of f^+l found hy equation (A3). 

The suggestion in reference 8 that the interval At he chosen s\affi- 
cently small so that there is little change in the values of tj.+l 

obtained from these two equations is followed in the present paper. 


Initial Values 

As tending of the column progresses, the strain in the convex flange 
eventually reaches a maximum and reverses. When reversal occurs, it is 
necessary to change the stress -strain relation for the convex flange 
which unloads elastically. Changing the stress -strain relation alters 
the form of g(f,r); th\is, the problem may he considered as consisting 
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of two parts, (l) prereversal and (2) postreversal, both of which require 
three initial valiies in order to start the recurrence process. 

Frereversal .- The required initial values for the prereversal phase 
can be determined by expanding the dliaensionless lateral deflection in a 
Taylor's series about t = 0, subject to the assumed initial condi- 
tions (1). The ratios of the moduli and Eg^^^E which appear 

in equation (21) are in5)licit functions of f and therefore differen- 
tiation of the equation is impractical .• However, during the initial, 
stage of loading the values of these ratios often may be assumed to be 
constant and eqiial to unity. Such an assumption permits differentiation 
of equation (21) in order that the higher -order derivatives involved in 
the Taylor's series may be found. The following relation obtained from 
the series expansion and the assimied Initial conditions provides the 
required initial values of f : 



efl2(2 + e2) 

2(^e/^i)^ 


51 


m 


For those stress -strain curves which have a small, or no, elastic 
region (for exanple, n = 2 in the Ramberg-Osgood stress -stra in curve 
representation), the assunptlon that the ratios of the moduli are unity 
is, of course, not valid. Nevertheless, the initial values of f 
obtained from equation (A^t-) may be xised eis the first trial values for 
an iteration procedure. An iteration expression can be derived by 
assuming that a second-degree parabola is sufficiently accurate to 
describe f" in the neighborhood of t = 0. Integration of the parab- 
ola yields the following expression for f in terms of f": 


f - 

- 21 ^ 


(n - 2)f"2 + 2{h - n)f"i + ^n - 6 + ^jf"^(Ar)^ + nf'^ At + f^ 

(n = 0, 1, 2) (A5) 


The trial values of f obtained from the series expansion are sub- 
stituted into equation (Al) to find the corresponding trial values 
of f", which are then substituted into equation (A5) to yield better 
values of f . The procedirre is repeated until the values of f converge. 

Postreversal . - As was previously mentioned, the integration process 
must be started anew at the point where strain reversal occurs. The 
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values of f, f, and f" at must be determined for they are 

used as Initial conditions for the new integration problem. Strain 
reversal will generally occur at some value of t which is not a mul- 
tiple of At, and, therefore, the point at which reversal occinrs must 
be located numerically or graphically. The required initial values 
of f and f" are found by the following procedure: 

1 . Find "^rev numerically or graphically from 

plots of against t. 

2 . Solve for f^.^^ by \ising equation (l6b). 

3. By using the known value of solve equation (Al) for f'^e- 

Tnfljmm ch as the velocity f ' of lateral deflection does not appear 
in the solution, it is more difficult to find i*'rev* ^ approximation 
to f 'rev obtained as follows: 

f'rev = f'a+/ 

J T=a 


where a is some val\je of t in the prereversal phase of the problem. 
By assuming that 


f ’ 


At 


(A7) 


equation (a 6) can be written as 


J p ~ * rsv 

f’dT (AB) 

T=a 

The first differences Af which are used in the solution are plotted 
at the midpoints of the intervals. The q-uantity (AT)f'j.ey is evaluated 
by numerically integrating the f " values between the last midinter- 
val point before reversal (t = a) and and adding the result to 

the value of Z 5 f at that last midinterval point. In the present case. 
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the error involved in the first-order approximation, equation (Aj), is 
small, and the use of the approximation is therefore justified. The 
errors in the term (Ar)f and the values of f calculated from 

equation (A5) are, of course, much smaller than the error in equation (A7). 

In addition to the initial conditions, two additional sets of f 
eind f" are required to restart the integration. These are obtained 
by using equation (A5) and the associated iteration procedure. The 
first trial values in this case are obtained by interpolating the pre- 
reversal solution, which was extended a few increments in t beyond 
the point at which reversal occurred. 


Inteiwal of Argument 

In the present paper the Interval At has been selected to pro- 
vide a m1.ni.Tmjm of 10 Intervals per period of the first natural fre- 
quency in transverse bending. For the faster loading rates, in which 
cases there is only a fraction of a period of oscillation before the 
maximum load is reached or before the strains become prohibitively 
large, as majny as 600 intervals per period of natural frequency were 
used. 
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TABLE 2.- BAWGE OF VELOCITIES EMPLOYED 


l^or velocity of compression-wave propagation 
= 16,7^ ft/sec and ^ = 56*^ 


a 

Velocity 

ft/sec 

in. /min 
(a) 

1 

8.14 

5861.0 

10 

2.57 

1855.0 

io2 

.814 

586.1 

io5 

.257 

185.5 

10^ 

.0814 

58.61 

lo5 

.0257 

18.55 

10^ 

.00814 

5.861 


^Maximum platen speed for commonly used 
hydraulic testing machines is approximately 
5 in. /min. 







■ -^ueB of eM displacement. 


f^ 




U(y^ ill vavH 
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Figure 4.- Variation of load with shortening for Initially curved Idealized 
H-section colunins with various rates of end displacenuent . n = 10. 
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Figure 5 *- Variation, of load with lateral deflection for dynamically 
loaded idealized E-aection columns vltti various Initial ctcrvatures 
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Figure 6.- Variation of Load vlth Bhortenlng for dynamically loaded 
idea3J.zed H-section columns vlth variaus initial curvB.ture8. 
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1 Figure 7-- Variation of load with Lateral deflection for dynamically 

i loaded, initially curved, Idealized H-eectlon col\imns of various 

lengths . 
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Figure 8.- Varlatioai of load vith lateral deflection for Initially curved 
idealized H- section columns vith various rates of end displacement, n = 2 
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Figure 9*- Variation of load witli shortening for initially curved idealized 
H- section column s 'with various rates of end displacement, n = 2. 
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